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ABSTRACT 

We investigate the properties of clusters of galaxies in the ACDM models with a step- 
like initial power spectrum. We examine the mass function, the peculiar velocities 
and the power spectrum of clusters in models with different values of the density pa- 
rameter Hq, the normalized Hubble constant h and the spectral parameter p, which 
describes the shape of the initial power spectrum. The results are compared with obser- 
vations. We also investigate the rms bulk velocity in the models, where the properties 
of clusters are consistent with the observed data. We find that the power spectrum 
of clusters is in good agreement with the observed power spectrum of the Abell-ACO 
clusters, if the spectral parameter p is in the range p — 0.6 — 0.8. The power spec- 
trum and the rms peculiar velocity of clusters are consistent with observations only if 
Qq < 0.4. The f2o = 0.3 models are consistent with the observed properties of clusters, 
if h = 0.50 - 0.63. For h = 0.65, we find that fl = 0.20 - 0.27. 
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1 INTRODUCTION 

The theory of galaxy formation based on gravitational in- 
stability describes how primordially generated fluctuations 
grow into galaxies and clusters of galaxies due to self-gravity 
of matter. The initial field of density fluctuations <5(x, t) can 
be decomposed into its Fourier components <5k(i) and ex- 
pressed in terms of the power spectrum P(k) — (|5k| 2 )- 

As the study of the large-scale structure in the Uni- 
verse has pushed to ever larger scales, several data sam- 
ples have suggested the presence of a peak in the power 
spectrum at the wavenumber k ~ 0.05ft Mpc -1 (or at the 
wavelength A ~ 120 - 130ft" 1 Mpc). Einasto et al. (1997) 
and Retzlaff et al. (1998) studied the spatial distribution of 
the Abell-ACO clusters and found that the power spectrum 
of the clusters has a well-defined peak at the wavenumber 
k = 0.052/i Mpc -1 . A similar peak in the one-dimensional 
power spectrum of a deep pencil-beam survey was detected 
by Broadhurst et al. (1990), and in the two-dimensional 
power spectrum of the Las Campanas redshift survey by 
Landy et al. (1996). The power spectrum of the spatial dis- 
tribution of APM galaxies also has a feature on the same 
scale (Caztanaga and Baugh 1998). Independent evidence 
for the presence of a preferred scale in the Universe at about 
130ft -1 Mpc comes from an analysis of high-redshift galax- 
ies. Broadhurst and Jaffe (1999) studied the distribution of 
the Lyman-break galaxies at redshift z ~ 3 and found a 5a 
excess of pairs separated by Az = 0.22 ± 0.02, equivalent 
to 130ft _1 Mpc for flat universe with the density parameter 
fi = 0.4 ±0.1. 



The power spectrum of density fluctuations depends on 
the physical processes in the early universe. The peak in the 
power spectrum of clusters at the wavenumber k ~ 0.05ft 
Mpc -1 may be generated during the era of radiation dom- 
ination or earlier. Standard cosmological models based on 
collisionless dark matter [e.g. cold dark matter (CDM)] and 
adiabatic fluctuations, when combined with power-law ini- 
tial power spectra, predict smooth power spectra of density 
fluctuations at z ~ 10 3 . The baryonic acoustic oscillations 
in adiabatic models may explain the observed power spec- 
trum only if currently favored determinations of cosmolog- 
ical parameters are in substantial error (e.g., if the density 
parameter f2o < 0.2ft; Eisenstein et al. 1998). 

One possible explanation for the observed power spec- 
trum of clusters is an inflationary model with a scalar field 
whose potential V(tp) has a local steplike feature in the first 
derivative. This feature can be produced by a fast phase 
transition in a physical field different from the inflaton field. 
An exact analytical expression for the scalar (density) per- 
turbations generated in this inflationary model was found by 
Starobinsky (1992). The initial power spectrum of density 
fluctuations in this model can be expressed as 



where the function S(k,p) can be written as 



(1) 



S(k,p) 



Hp - 1) 



fi(y) sin 2y + -cos2i/ 
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Figure 1. The function S(k,p)/p 2 for different values of the parameter p. For the models with p < 1 and p > 1, the step parameter was 
chosen to be ko = O.OWh Mpc -1 and fco = 0.03fo Mpc -1 , respectively. 



9(p-l) 2 / 2 (y) 
2y 2 



h (y) + fi (y) cos 2y - - sin 2y 



(2) 



Here, the function y — k/ko, fi(y) = 1 — y~ and f2{y) = 
l+y -2 . The initial power spectrum in this model depends on 
two parameters fco and p. The parameter fco determines the 
location of the step and the parameter p - the shape of the 
initial spectrum. For p = 1, we recover the scale-invariant 
Harrison-Zel'dovich spectrum (S(k, 1) = 1). At present, the 
initial spectrum (1,2) is probably the only example of an 
initial power spectrum with the desired properties, for which 
a closed analytical form exists. 

Fig. 1 shows the function S(k,p)/p 2 for different values 
of the parameter p. For the models with p < 1 and p > 1, the 
step parameter was chosen to be fco = O.OI6/1 Mpc -1 and 
fco = 0.03ft Mpc -1 , respectively. In this case, in the models 
with p < 1, the power spectrum has a well-defined maxi- 
mum at the wavenumber fc ~ 0.05ft Mpc -1 and a second 
maximum at fc ~ 0.1ft Mpc -1 . In the models with p > 1, 
the picture is inverted. The power spectrum has a flat up- 
per plateau at the wavenumbers fc < 0.05ft Mpc -1 , a sharp 
decrease on smaller scales (fc = 0.05 — 0.1ft Mpc -1 ) and a 
secondary maximum at fc ~ 0.15ft Mpc -1 . 

Lesgourgues, Polarski & Starobinsky (1998, hereafter 
LPS) compared the CDM models with a steplike initial spec- 
trum (1,2) with observational data. They studied the rms 



mass fluctuation on an 8ft -1 Mpc scale, as; the rms bulk 
velocity of galaxies and the cosmic microwave background 
(CMB) anisotropics on different angular scales. In this pa- 
per we continue this research and investigate the properties 
of galaxy clusters in these models. We examine the mass 
function, the peculiar velocities and the power spectrum of 
clusters for different values of the density parameter Qo , the 
normalized Hubble constant ft and the spectral parameter p. 
The results are compared with observations. We also inves- 
tigate the rms bulk velocity and as in the models, where the 
mass function, the peculiar velocities and the power spec- 
trum of clusters are consistent with the observed data. 

In their study LPS assumed that the parameter as lies 
in the interval a s = (0.57 ± O.OG)^ 56 . This interval was 
derived by White, Efstathiou & Frenk (1993) by analysing 
the mass function of clusters. For Qq — 0.3 this gives as = 
1.0 — 1.25. In this paper we examine the observed values 
of the mass function of galaxy clusters in more detail and 
obtain lower values for the parameter as- For Qq = 0.3 we 
find that as = 0.8 — 1.05. Therefore, the allowed values for 
the parameter p in the (fio, ft) plane that we obtain are 
different from those found by LPS. 

We examine flat cosmological models with the density 
parameter fi = 0.2 — 0.5 and the normalized Hubble con- 
stant ft = 0.5 — 0.8. These parameters are in agreement with 
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measurements of the density parameter (e.g. Bahcall et al. 
1999) and with measurements of the Hubble constant using 
various various distance indicators (e.g. Tammann 1998). 
To restore the spatial flatness in the low-density models, 
we assume a contribution from a cosmological constant: 
SIa = 1 — fio- The Hubble constant is written as Ho = 100ft 
km s -1 Mpc -1 . The transfer function T(k) is computed with 
the fast Boltzmann code CMBFAST developed by Seljak & 
Zaldarriaga (1996). The models are normalized using the 
COBE normalization derived by Bunn & White (1997). We 
assume that the initial density fluctuation field in the Uni- 
verse is a Gaussian field. In this case, the power spectrum 
provides a complete statistical description of the field. 

This paper is organized as follows. In Section 2 we study 
the mass function of clusters of galaxies and compare the re- 
sults with observations. In Section 3 we examine the peculiar 
velocities of galaxy clusters. In Section 4 we investigate the 
power spectrum of clusters. Discussion and summary are 
presented in Section 5. 



2 THE MASS FUNCTION OF CLUSTERS OF 
GALAXIES 

To study the mass function of clusters we use the Press- 
Schechter (1974, PS) approximation. The PS mass function 
has been compared with N-body simulations (Efstathiou et 
al. 1988; White, Efstathiou & Frenk 1993; Lacey & Cole 
1994; Eke, Cole & Frenk 1996; Borgani et al. 1997a) and 
has been shown to provide an accurate description of the 
abundance of virialized cluster-size halos. In the PS approx- 
imation the number density of clusters with the mass be- 
tween M and M + dM is given by 



mass Mi. 5, is related to the window radius R as 



n(M)dM 



Pb 



dcr(M) 



nMa 2 (M) dM 



8 2 t 



2a 2 (M) 



cxp - n _,,i^ dM. 

(3) 

Here pb is the mean background density and St is the linear 
theory overdensity for a uniform spherical fluctuation which 
is now collapsing; St = 1.686 for Qo = 1, with a weak de- 
pendence on Qo for flat and open models (Eke et al. 1996; 
Kitayama & Suto 1996). In this paper we will use values of 
St found by Eke et al. (1996) for flat models. The function 
a(M) is the rms linear density fluctuation at the mass scale 
M. We will use the top-hat window function to describe ha- 
los. For the top-hat window, the mass M is related to the 
window radius R as M = 4irpbR [i /3. In this case, the number 
density of clusters of mass larger than M can be expressed 



n c i(> M) 



-f 

J M 



n{M')dM' = 



(2tt)3/2 



f-OO 

Jr a 



St da(r) 
2 (r) dr 



exp 



2a 2 (r) 



Fig. 2 shows the cluster mass function in the flat mod- 
els with fi = 0.3 and h — 0.65. The threshold density 
8 t = 1.675. We investigated the cluster masses within a 
1.5ft -1 Mpc radius sphere around the cluster center. This 



R = 8.43!X; 



Mi., 



6.99 x 10 14 a>ft -1 M G 



(ft -1 Mpc). 



(5) 

Here the parameter a describes the cluster mass profile, 
M(r) ~ r a , at radii r ~ 1.5ft -1 Mpc. Numerical simula- 
tions and observations of clusters indicate that the param- 
eter a 0.6 — 0.7 for most of clusters (Navarro, Frenk & 
White 1995; Carlberg, Yee & Ellingson 1997). In this paper 
we use a value a = 0.65. 

Fig. 2 shows also the mass function of clusters of galax- 
ies derived by Bahcall and Cen (1993, BC) and by Girardi 
et al. (1998, G98). BC used both optical and X-ray observed 
properties of clusters to determine the mass function of clus- 
ters. The function was extended towards the faint end using 
small groups of galaxies. G98 determined the mass function 
of clusters by using virial mass estimates for 152 nearby 
Abell-ACO clusters including the ENACS data (Katgert et 
al. 1998). The mass function derived by G98 is somewhat 
larger than the mass function derived by BC, the difference 
being larger at larger masses (see Fig. 2). 

Let us consider the amplitude of the mass function of 
galaxy clusters at Mi. 5 = 4 • 10 14 ft -1 M Q . For this mass, the 
cluster abundances derived by BC and G98 are n(> M) = 
(2.0 ± 1.1) • 10~ 6 ft 3 Mpc" 3 and n(> M) = (6.3 ± 1.2) • 
10 -6 ft 3 Mpc -3 , respectively. By analysing X-ray properties 
of clusters, White, Efstathiou & Frenk (1993) found that 
the number density of clusters with the mass Mi. 5 « 4.2 • 
10 14 ft -1 M e is n(> M) = 4 • 10~ 6 ft 3 Mpc~ 3 . 

We derived the limits for the parameter p, assum- 
ing that the mass function of galaxy clusters at Mi. 5 = 
4 • 10 14 ft -1 M o is in the range (2 - 6.5) • 10~ 6 ft 3 Mpc" 3 . 
Fig. 3a shows the results for the models with ft = 0.65 for 
different flo and Fig. 3b for the models with f2o = 0.3 for 
different ft. For the model with Qo = 0.3 and ft — 0.65, 
we find that p — 0.79 — 1.0. Fig. 2 demonstrates the mass 
function of clusters in this model for p — 0.79 and p = 1.0. 

LPS used high values of the parameter as and found 
that one of the best-fit models is the model with parame- 
ters Oo = 0.3, ft — 0.7, p — 0.8. We find that the number 
density of clusters in this model is substantially higher than 
observed; for the mass Mi. 5 = 4 • lO 14 ft _1 M , the cluster 
abundance n(> M) = 9.8 x 10 _6 ft 3 Mpc" 3 . For fi = 0.3 
and h = 0.7, the mass function of clusters is consistent with 
the observed data, if p = 0.88 — 1.13. 



3 PECULIAR VELOCITIES OF CLUSTERS OF 
GALAXIES 

The observed rms peculiar velocity of galaxy clusters has 
been studied in several papers (e.g. Bahcall, Gramann & Cen 
1994, Bahcall and Oh 1996, Borgani et al. 1997b, Watkins 
1997). In this paper we use the results obtained by Watkins 
(1997). He developed a likelihood method for estimating the 
rms peculiar velocity of clusters from line-of-sight velocity 
measurements and their associated errors. This method was 
applied to two observed samples of cluster peculiar veloc- 
ities: a sample known as the SCI sample (Giovanelli et al. 
1997) and a subsample of the Mark III catalog (Willick et al. 
1997). Watkins (1997) found that the rms one-dimensional 
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Figure 2. The cluster mass function in the models with Q.q = 0.3, h = 0.65, p = 1.0 (solid line) and p = 0.79 (dashed line). Open circles 
and squares show the mass function of galaxy clusters derived by Bahcall and Cen (1993) and by Girardi ct al. (1998), respectively. The 
open triangle is the result obtained by White, Efstathiou & Frenk (1993). 



cluster peculiar velocity is 265l7g 6 km s 1 , which corre- 
sponds to the three-dimensional rms velocity 459ti|o km 

To investigate the peculiar velocities of clusters in our 
models, we use the linear theory predictions for peculiar 
velocities of peaks in the Gaussian field. The linear rms ve- 
locity fluctuation on a given scale R at the present epoch 
can be expressed as 

a v (R) = H f{n )a- 1 (R), (6) 

where f(Oo) ~ fif!' 56 is the linear velocity growth factor in 
the flat models and Oj is defined for any integer j by 

a) = ^_ J P(k)W'\kR)k 23+2 dk. (7) 

Bardeen et al. (1986) showed that the rms peculiar velocity 
at peaks of the smoothed density field differs systematically 
from a v (R), and can be expressed as 

a p (R) = a v (R) y/l - a^/afal ~ (8) 

Suhhonenko & Gramann (1999) examined the linear 
theory predictions for the peculiar velocities of peaks and 
compared these to the peculiar velocities of clusters in N- 
body simulations. The N-body clusters were determined as 
peaks of the density field smoothed on the scale R ~ 1.5/i _1 



Mpc. The numerical results showed that the rms peculiar ve- 
locity of small clusters is similar to the linear theory expecta- 
tions, while the rms peculiar velocity of rich clusters is higher 
than that predicted in the linear theory. The rms peculiar ve- 
locity of clusters with a mean cluster separation d cl = 30ft _1 
Mpc was ~ 18 per cent higher than that predicted by the 
linear theory. We assume that the observed cluster sample 
studied by Watkins (1997) corresponds to the model clus- 
ters with a separation d c i ~ 30/i -1 Mpc (n ci ~ 3.7 • 10 _5 /i 3 
Mpc -3 ) and determine the rms peculiar velocity of the clus- 
ters, v c i, as 

Vd = 1.18 a p ( R) , (9) 

where the radius R — 1.5ft _1 Mpc. 

Fig. 3 shows the limits for p in different models obtained 
on the basis of Watkins (1997) results. Fig. 3a shows the 
results in the h = 0.65 models for different fio and Fig. 3b 
in the Qo = 0.3 models for different h. In the region studied 
in Fig. 3b, the rms peculiar velocity is larger than 329 km 
s _1 . For the model with fio = 0.3 and h = 0.65, we find that 
the initial parameter p > 0.87. 

For high values of fi and h, peculiar velocities are 
larger than observed for any values of the parameter p. The 
velocities are sensitive to the amplitude of the large-scale 
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Figure 3. The limits for p in the ACDM models with a stcplikc initial power spectrum. Solid lines show the constraints obtained from 
the mass function of clusters and dashed lines show the constraints obtained by analysing the peculiar velocities of clusters, (a) h = 0.65. 
(b) Qq = 0-3. In the region studied peculiar velocities are larger than 329 km s _1 . 



fluctuations at wavenumbers k < O.lh Mpc -1 . By increas- 
ing p, if p > 1, the power spectrum on these wavenumbers 
and, therefore, the velocities remain almost unchanged (see 
Fig. 1). Fig. 3a shows that in the h — 0.65 models with 
Qa > 0.35, the peculiar velocities are not consistent with 
observations. For the flo = 0.4 model with h = 0.5 and 
h = 0.6, we find that peculiar velocities are consistent with 
observations, if p > 0.86 and p > 1.26, respectively. 

If we compare the observational constraints obtained by 
studying the mass function and peculiar velocities of clus- 
ters of galaxies, we see from Fig. 3 that the mass function 
and the peculiar velocities of clusters are consistent with the 



observed data only in a small interval of the parameter p. 
For the model with flo — 0.3 and h — 0.65, we find that 
p = 0.87 - 1.0. For the fi = 0.3 models with h = 0.5 and 
h = 0.6, we find that p = 0.57 - 0.64 and p = 0.76 - 0.88, 
respectively. 



4 THE POWER SPECTRUM OF CLUSTERS 
OF GALAXIES 

Let us now consider the power spectrum of clusters in these 
models where both the mass function and the peculiar ve- 
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b ci = i- 



locities are consistent with observations. To investigate the 
power spectrum of clusters, we can also use the PS approach. 
The power spectrum of clusters for a given number density 
in the PS approximation can be expressed as (Gramann & 
Suhhonenko 1999, hereafter GS) 

Pd(k) = b 2 cl P(k), (9) 

where the cluster bias parameter b c i is 

3 Z* 00 1 da(r) 2 2 dr 

(10) 

Here, the function y = 8 t /a(r) and a(r) is the rms linear 
density fluctuation at the radius r. The cluster bias param- 
eter bd depends on the minimal mass M (or the window 
radius 7?) of clusters and on the power spectrum of density 
fluctuations, P(k), which determines the function a(r). For 
fixed P(k) and n c i, the minimal mass M (or scale R) can be 
determined by inverting equation (4) . 

Observations provide the distribution of clusters in red- 
shift space, which is distorted due to peculiar velocities of 
clusters. On large scales, where linear theory applies, the 
power spectrum of matter density fluctuations in redshift 
space is given by (Kaiser 1987): 



P°(k) = 



1 | 2/(n„) | Ago) 



p(k). 



(ii) 



Using the PS approximation (10), relation (11) takes the 
form 



pm = 



1 | 2/(n ) | / 2 (n ) 



36 c 



5% 



KiP(k). (12) 



Equation (12) determines the power spectrum of clusters for 
a given n c i in redshift space. 

GS examined the power spectrum of clusters in the 
PS theory and in N-body simulations. They determined 
the power spectrum of clusters for mean separations d c i = 
30 — 40ft -1 Mpc. The numerical results showed that at 
wavenumbers k < O.l/i -1 Mpc, the power spectrum of clus- 
ters in the simulations is linearly enhanced with respect to 
the power spectrum of the matter distribution. However, the 
amplitude of the spectrum of clusters was somewhat lower 
than predicted by the approximation (12). It is possible that 
the linear approximation (12) overestimates the power spec- 
trum of clusters due to dynamical effects that are not taken 
into account in this approximation. The power spectrum of 
clusters can be expressed as 



pm = f 



1 | 2/(Oo) | f{n ) 



Sb c , 



b 2 cl P(k), (13) 



where the factor F = 0.7—0.8. The factor F depends slightly 
on the model. GS examined the distribution of clusters 
in two cosmological models which start from the observed 
power spectra of the distribution of galaxies and clusters of 
galaxies. In the model (1), the initial linear power spectrum 
of density fluctuations was chosen in the form P(k) oc k~ 2 
at wavelengths A < 120fo -1 Mpc. In the model (2), GS as- 
sumed that the initial power spectrum contains a primordial 
feature at the wavelengths A ~ 30 — 60h -1 Mpc. They found 
that F = 0.8 and F = 0.7 in the model (1) and model (2), 
respectively. In this paper we use a value F = 0.75. 

Fig. 4 shows the redshift-space power spectrum of the 



model clusters with a mean separation d c i = 34h~ Mpc. For 
comparison, we show the power spectrum of the Abell-ACO 
clusters determined by Einasto et al. (1999). This spectrum 
represents the weighted mean of the power spectra deter- 
mined by Einasto et al. (1997) and Retzlaff et al. (1998). 
Einasto et al. (1997) determined the power spectrum of the 
Abell-ACO clusters from the correlation function of clusters, 
while Retzlaff et al. (1998) estimated the power spectrum di- 
rectly (see Einasto et al. 1999 for details). The power spec- 
trum of the distribution of the Abell-ACO clusters peaks at 
the wavenumber k = 0.052/i Mpc -1 . The mean intercluster 
separation of the Abell-ACO clusters is d ct ~ 34/?, -1 Mpc 
(n r j ~ 2.5 ■ 10~ 5 fr 3 Mpc -3 ) (Einasto et al. 1997, Retzlaff et 
al. 1998). 

Fig. 4 shows the power spectrum of clusters predicted 
in the H = 0.3 models with (h = 0.5, p = 0.6), (h = 0.6, 
p = 0.8), (h = 0.65, p = 0.95) and (h = 0.7, p = 1.05). For 
these models, the mass function and the peculiar velocities 
are consistent with the observed data (see Fig. 3b). The 
models with p — 0.6 and p — 0.8 are in good agreement 
with the observed power spectrum, while the models with 
p — 0.95 and p = 1.05 are not consistent with the observed 
data. In the models with p > 0.95, we do not see a peak in 
the power spectrum at the wavenumber k ~ 0.05/i Mpc -1 . 

We used also the \ 2 test to calculate the probability 
that the models fit the observed power spectrum. The ob- 
served power spectrum represents the weighted mean of the 
power spectra determined by Einasto et al. (1997) and Ret- 
zlaff et al. (1998) and it is not clear how many points of the 
power spectrum are independent. As a first step we used 
all the data points in Fig. 4. The models with p = 0.6 and 
p — 0.8 are consistent with the observed power spectrum at 
a confindence level higher than 90%. For the models with 
p = 0.95 and p = 1.05, the probability that the models fit 
the observed power spectrum is less than 30%. 

We studied also the power spectrum of clusters in other 
models with different values of (fio, h) and found similar re- 
sults. The power spectrum of clusters is in good agreement 
with the observed power spectrum of the Abell-ACO clus- 
ters, if the initial parameter p is in the range p — 0.6 — 0.8. In 
the models with h — 0.65, the allowed values for the param- 
eter p are in this range, if fio = 0.20 — 0.27 (see Fig. 3a). Sim- 
ilarly, the Qq — 0.3 models are consistent with the observed 
mass function, peculiar velocities and the power spectrum 
of clusters, if h = 0.50 - 0.63 (Fig. 3b). 

In the fio = 0.4 model with h — 0.5, the peculiar veloci- 
ties are larger than observed, if p < 0.86. For higher values of 
fio and h, this limit for p increases. Therefore, for fio > 0.4 
and h > 0.5, the observed power spectrum and peculiar ve- 
locities of clusters are not consistent with each other. Either 
the observed peculiar velocities are underestimated, or the 
observed peak in the power spectrum of clusters is overes- 
timated. The peculiar velocities and the power spectrum of 
clusters are consistent with observations only if fi < 0.4. 



5 DISCUSSION AND SUMMARY 

In this paper, we have examined the properties of clusters of 
galaxies in the ACDM models with a steplike initial power 
spectrum (1,2) that depends on two parameters ko and p. 
The parameter ko determines the location of the step and 
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Figure 4. The power spectrum of the model clusters and the Abell-ACO clusters (open circles) with a separation d c i = 34h~ 1 Mpc. 
The heavy lines show the power spectra of clusters determined by approximation (13) and the light lines the corresponding linear power 
spectra of matter fluctuations. We have shown the spectra for the S7o = 0.3 models with h = 0.5, p = 0.6 (dot-dashed lines), h = 0.6, 
p = 0.8 (solid lines), h = 0.65, p = 0.95 (dotted lines) and h = 0.7, p = 1.05 (dashed lines). 



the parameter p - the shape of the initial spectrum. For the 
models with p < 1 and p > 1, the step parameter was chosen 
to be ko = 0.016h Mpc -1 and fco = 0.03/i Mpc -1 , respec- 
tively. We investigated the mass function, peculiar velocities 
and the power spectrum of clusters in models with different 
values of f2o, h and p. 

We found that the mass function and the peculiar ve- 
locities of clusters are consistent with the observed data only 
in a small interval of the parameter p. For the model with 
Qo = 0.3 and h = 0.65, we find that p = 0.87 - 1.0. For the 
Qo = 0.3 models with h — 0.5 and h = 0.6, we find that 
p — 0.57 — 0.64 and p = 0.76 — 0.88, respectively. 

The power spectrum of clusters in the ACDM models 
with a steplike initial power spectrum is in good agreement 
with the observed power spectrum of the Abell-ACO clus- 
ters, if the initial parameter p is in the range p = 0.6 — 0.8. 
In the models with h = 0.65, the allowed values for the 
parameter p are in this range if flo = 0.20 — 0.27. The 
0,0 = 0.3 models are consistent with the observed mass func- 
tion, peculiar velocities and the power spectrum of clusters, 
if h = 0.50 — 0.63. The peculiar velocities and the power 
spectrum of clusters are consistent with observations only if 
O < 0.4. 

We also studied the rms bulk velocity in the models, 



where the mass function, the peculiar velocities and the 
power spectrum of clusters are consistent with observations. 
We choosed three models in the allowed parameter space 
and studied these models in more detail. Tabel 1 lists the 
parameters of the models studied. We have determined the 
number density of clusters with mass Mi. 5 = 4- 1O 14 /i _1 M0, 
n c i \ the rms mass fluctuation on an 8h _1 Mpc scale, as; the 
rms peculiar velocity of clusters, v c i , and the rms bulk veloc- 
ity for a radius r — 50h~ Mpc, V50. The rms bulk velocity 
was determined by using eq. (6) . The observed bulk veloci- 
ties are determined in a sphere centered on the Local Group 
and represent a single measurement of the bulk flow on large 
scales. The observed bulk velocity derived from the Mark III 
catalog of peculiar velocities for r — 50ft" 1 Mpc is 375 ± 135 
km s _1 (Kolatt & Dekel 1997). In the models studied, the 
rms bulk velocity is ~ 270 — 285 km s _1 , which is consistent 
with the observed data. 



Therefore, in many aspects the ACDM models with a 
steplike initial spectrum fit the observed data. Further work 
is needed to study the properties of galaxies and clusters of 
galaxies in these models in more detail. 
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Table 1. Results of the tests for the different models. 

Q h p n c i a s v ct Vso 

(l(r 6 h 3 Mpc~ 3 ) (km s- 1 ) (km s" 1 ) 

0.20 0.75 0.70 2.9 1.05 610 270 

0.25 0.65 0.75 2.8 0.95 610 270 

0.30 0.60 0.80 3.3 0.90 635 285 
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